Directed coherent transport due to Bloch oscillation in two dimensions 
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We point out that in higher dimensions, in contrast to the one dimensional case considered usually, 
Bloch oscillation driven by a static force can induce transport of the wave packet. The wave packet 
oscillates constantly, but on a larger time scale it drifts at a constant velocity permanently. As a 
noteworthy feature, the net transport in the long run is always normal to the external force and 
thus controlled by it. We verify this prediction numerically and discuss its experimental realization 
both with cold atoms in optical lattices and with two dimensional photonic lattices. 
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Bloch oscillation is a peculiar response of a particle in 
a periodic potential to an external force Under a 
static uniform force, the particle perforins an oscillatory 
motion in the real space without falling down to the po- 
tential minimum at infinity as in free space. This counter- 
intuitive phenomenon is purely quantum and has its root 
deep in the band structure associated with the periodic 
potential. The periodic potential mixes the plane waves 
into Bloch waves, which are classified by band number n 
and wave vector fc. For a weak external force, the gap 
between the bands protects the particle from transition 
into other bands, with the only effect being that the wave 
vector is dragged across the Brillouin zone (BZ). The os- 
cillation of the center-of-mass of the wave packet is then 
related to the fact that Bloch waves are periodic with 
respect to fc. 

So far, Bloch oscillation has been investigated and ob- 
served in a variety of systems, such as semiconductor 
superlattice [3] , cold atoms in optical lattices [1, 01 , and 
photonic lattices [B-0|- However, most of these works are 
confined to one dimension (some exceptions being [1, 0]). 
In one dimension, the Brillouin zone [—G/2, +G/2] has 
the topology of a circle since — G/2 is identified with 
+G/2. Therefore, in the fc-space, the motion of the par- 
ticle is simple — it transverses the Brillouin zone repeat- 
edly. Quantitatively, we have the Newton second law 
type equation of motion Q : 



dk ^, , 



(1) 



where F{t) is the time-dependent external force, and the 
wave vector k is identified with k -\- nG for an arbitrary 
integer n {h — 1 in this paper). In the static case F{t) = 
F, after a period of T = G/F, the wave vector returns to 



its initial value ko and so does the center-of-mass [sl, Q : 
dr 

F / dt— = F dtVkE{ko + Ft) 

Jo "'0 

= E{ko + G) - E{ko) = 0. (2) 

Here E{k) is the dispersion function of the relevant en- 
ergy band. Therefore, in one dimension, Bloch oscilla- 
tion induced by a static force does not lead to trans- 
port. To induce a net transport, a time-dependent force 
is needed. This makes sure that the Brillouin zone is 
transversed non-uniformly and thus the center-of-mass 
displacement accumulated in different /c-regions do not 
cancel each other. This is essentially the scheme em- 
ployed to induce the so-called super Bloch oscillations 
and macroscopic transport with cold atoms in one di- 
mensional optical lattices [lol - [l2l |. 

In this paper, we would like to point out that the situa- 
tion is different in higher dimensions. In a higher dimen- 
sional potential, suppose the wave vector returns ever 
to its initial value by been dragged by some reciprocal 
lattice vector G, the counterpart (or generalization) of 
Eq. © is 



, dr 
dt — 

dt 



[ d{Ft)-\/kE{ko + Ft) 
Jo 

Eiko + G)- E{ka) = 0. (3) 



Thus we can only arrive at the conclusion that as the 
wave vector returns to its initial value, the displacement 
of the center-of-mass of the wave packet must be perpen- 
dicular to the external force. But it does not necessarily 
vanish. 

In the following, we give an explicit expression of the 
displacement in one period. We take the tight binding 
model to consider the Bloch oscillation. Note that, in 
a real lattice, Bloch oscillation is a good picture only 
when the single-band approximation is valid, i.e., when 
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the band-band transition or the Zener tunnelhng is neg- 
hgible. With a tight binding model on a simple Bravais 
lattice, there is one and only one band, and the system 
is then free of band-band transition. Conversely, if we 
restrict to a single band, then effectively we are dealing 
with a tight binding model on a simple lattice. The Wan- 
nier functions associated with that band in each primitive 
cell serve as the basis states on each site. Of course, by 
restricting to a single band, we necessarily misses the po- 
tential Berry phase effect and the transport associated 
p^ . However, since generally the Berry curvature is 
small in most area of the BZ [IJ], or even vanishes iden- 
tically in the presence of both time reversal and spatial 
inversion symmetry [isl. [l5l|. we neglect it in this paper. 

One advantage of shaping the problem on a lattice is 
that, quantities of only a metric value are irrelevant. To 
be precise, the exact shape of the lattice, e.g. the lengths 
of the two basis vectors and the angle between them are 
of no concern. We then simply take the two basis vectors 
along the x and y directions, and assume their lengths to 
be unity. The first Brillouin zone is then [— tt, tt] x [— tt, tt]. 
The original real lattice differs from this square lattice by 
only an affine transform. 

The tight binding Hamiltonian is Htb — 
m Jrn\l){l + m]- Here I and m both take values 
among all the integer pairs (7711,7712), ~oo < mi^2 < 00, 
and \l) denotes the Wannier funtion at site I. Because 
Htb = H^g, we have Jm — J-m- The Bloch states 



are \E'fc ~ X^m^ '"l"^)' 'where k = (fci,fc2) is the wave 
vector. They are eigenstates of Htb with eigenvalues 
E{k) = — J2m "^rn Gxp{ik ■ m) . In many cases, the 
original Hamiltonian is time reversal invariant, and 
therefore E{k) — E{—k) necessarily. This implies 
Jm = J-m- Combined with the previous condition, this 
implies is real. 

Now suppose a linear potential is applied to the system. 
The Hamiltonian is now H — Htb — J2m ^ ■ rnlm) (m|. 
In the semi-classical theory, the details of the wave packet 
are neglected and it is characterized just by a pair of 
variables (t, fe), which are respectively the average values 
of the position and moment of the wave packet. The 
semi-classical equations of motion are J3\ 



dk 

dr 



F, (4) 
VE{k) ~ —i}^ mJ^ exp(ifc ■ m). (5) 



These equations can be understood as follows. First note 
that an initial Bloch state remains as a Bloch state all 
the time, more precisely, the time varying Bloch state 
\l/fe(,+Fte~'-/^3 solves the Schrodinger equa- 
tion idijj/dt ~ Hip. This actually implies or solves 
Eq. Now for an initial wave function V'm(O) = 

/^^ dkf{k)e''' '^, where /(fc) = /(fco + Sk) = g{6k) is a 
function narrowly localized around fco (so that it is mean- 
ingful to say the wave vector of the initial wave function 



is fco) 5 the wave packet at time t is 
Jbz 

' d{5k)g{5k)e'^''<^-io dt,VE(k,+Ft,)) 

^^i{ko+Ft)-m^-i f^'^ dtiE{ko+Fti) ^ ^g^ 

Here the phase factor in the first line is expanded around 
fco. This is legitimate because the main contribution to 
the integral comes from a small neighborhood of feg . The 
second line implies that the wave packet is translated 
forward by 



r{t) - r(0) = [ dtiVE{ko + Fh) 
Jq 

= V mJ„ / dtie'^'"' 
Jo 



m+F-mt 



which solves Eq. ([S]). 

Assume that there exist some time T such that FT is 
equal to a reciprocal lattice vector, that is. 



FT = {FiT, F2T) = (27r(7, 27rr), 



(8) 



where q and r are some co-prime integers. For F1/F2 = 
q/r, the explicit value of T is T = 27rq/Fi. After one 
period of T, the wave vector returns to its initial value 
since fc is identified with k+{27rq, lixr). The displacement 
of the wave packet in this period is calculated as 

/■T 

= mJ^sin(m • fco)(5F m- (9) 

^ — ^m 

This equation is our central result. The Dirac function 
means that for a given F in some direction, only those 
m which are perpendicular to F contribute to the net 
displacement. Physically, this means in the long run the 
wave packet drifts perpendicular to the external force. 
This is quite reasonable. If after one period of T, the 
wave packet is shifted somewhat along the direction of 
the external force, then in the long run, the shift of the 
wave packet along the direction of the external force will 
diverge to A-co or —00 linearly with time. As the system 
is conservative, the kinetic energy of the particle will go 
to -l-oo or —00 respectively. However, this is impossible 
since in a definite band, Eik) has both upper and lower 
bounds. The argument applies equally well to any dimen- 
sions and holds as long as the single band approximation 
is valid. 

Another remark is worthy. Our assumption ^ about 
the direction of F seems quite stringent. Indeed, since 
the measure of rational numbers is zero, for a generic F, 
the ratio F1/F2 is irrational and thus Eqs. ^ and © 
seem irrelevant. Moreover, even restricted to the field 
of rational numbers, T and Dt are pathologically sensi- 
tive to the ratio F1/F2 (e.g., suppose F1/F2 is perturbed 
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from 1/2 to 101/200, T enlarges by a factor of 100 while 
Dt shrinks almost to zero since the associated hopping 
coefficients are exponentially small) . The way out of this 
anomaly is that T and Dt both refer to the overall be- 
havior in the long run and thus they become irrelevant if 
T ^ T, the observation time span. We further note that 
in a given finite time interval, the displacement r{t)—r{0) 
in Eq. ([7]) is a well behaved smooth function of F1/F2. 
Thus, as long as the motion in a given time interval is 
concerned, the concerns above are irrelevant. Of course, 
to make the directed transport effect as significant as pos- 
sible, it is preferable to choose F in the vicinity of those 
directions (q, r) with q and r not large. 

In the hindsight, we see from Eq. ^ why Bloch oscil- 
lations in one dimension and in higher dimensions behave 
differently. In one dimension, all k vectors undergo the 
same trajectory and thus the phases accumulated are the 
same for all the k vectors. This means after one period, 
the wave function itself returns to its initial value (up to 
a global phase), not only its center-of-mass. In contrast, 
in higher dimensions, different k vectors may take differ- 
ent (parallel) trajectories in the fe-space, though they all 
return to their initial values simultaneously. The phases 
accumulated in one period for k vectors on different tra- 
jectories differ in general and this dispersion leads to the 
shift of the center-of-mass of the wave packet. 

In the following we apply the formalism developed 
above to some concrete case. We consider the scenario of 
cold atoms in a planar triangular optical lattice. This lat- 
tice can be constructed by interfering three laser beams 
travelling 27r/3 with respect to each other in the xy plane 
and all polarized in the z direction fl^ . The laser beams 
are all of wave length A and amplitude E. Suppose the 
laser is red detuned from the atomic transition, then the 
minima (of value Vq oc 9E^) of the optical potential form 
a triangular lattice (see Fig.[lja)). By solving the energy 
bands exact numerically and then performing Fourier 
transform by using the expression of E(k), we can solve 
the hopping amplitudes Jm which are controlled by Vq. 
In Fig. [TJb) , we show the logarithm of their magnitudes 
In I J^l in the case of Vq = —l.5Er, with Er = /2m\^ 
being the recoil energy. We note that the site-site hop- 
ping decreases monotonically with the site-site distance. 
Moreover, the Dq point group symmetry of the lattice is 
respected perfectly. In our numerical simulations below, 
we will only preserve site-site hoppings belonging to the 
first three hopping values Ji (i = 1,2,3) according to 
their magnitudes. In the specific case of Vq = —\.hEr, 
(Ji, J2, J3) = (0.0765, -0.0149, -0.0078)£;r. Those site- 
site hoppings belonging to Ji are J±(i.o), >^±(o,i): and 
J_j-(i i), and those belonging to J2 are J±(2.i)i •^±(1.2)7 
and J±(_i,i), and to J3 are J±(2,o), -^±(0,2), and J±(2,2)- 

The initial wave packet is assumed to be a Gaussian, 
i'mit = 0) = Aexp(-(rn^ _|_ m2)/(T^ + iklmi + iklm2), 
with m — (mi,m2). Here a = 20 is the width of the 
wave packet and fco — {kg, kg) = (0.05, 0.03) is the initial 
wave vector, while A is some normalization factor. The 
width cr 3> 1 so that the wave function is well localized 



around feg in the fc-space. Three cases of different forces 
are investigated. They are (i) F/Ji — (0.5,-0.5), (ii) 
F/Ji = (0.7,-0.7), and (iii) F/Ji = (0.4,-0.8), respec- 
tively. The Schrodinger equation id')jj /dt = is solved 
on a 121 X 121 lattice by the fourth order Runge-Kutta 
method. The time evolutions of the center-of-mass of 
the wave packet (mi) and (7712), are shown in Fig. [2ja). 
As another respect, the trajectories (mi) versus (m2) are 
shown in Fig. ^h) . The semi-classical predictions accord- 
ing to Eqs. (|H) and (O agree with the exact numerical 
calculations perfectly and thus are not shown. In all the 
cases, we see that though there are temporary oscilla- 
tory deviations, in the long run, the wave packet drifts 
at a constant velocity perpendicular to the external force. 
In cases (i) and (ii), the drifting direction is along (1,1), 
while in case (iii), it is along (2, 1), as evident in Fig.[2Ub). 
We also note, by comparing case (i) with case (ii), that 
the oscillation amplitude is suppressed for a larger force, 
and on the other hand, the oscillation period is short- 
ened. These features are reminiscent of Bloch oscillation 
in one dimension. For a larger force, the Brillouin zone is 
transversed faster and thus the period is shortened, and 
so is the oscillation amplitude because the deviation ac- 
cumulated decreases accordingly. The three cases show 
how the transport direction and the detailed temporary 
motion of the wave packet can be controlled by the di- 
rection and magnitude of the external force, respectively. 

In conclusion, we have demonstrated that Bloch oscil- 
lations in higher dimensions, in contrast to their coun- 
terpart in one dimension, can lead to transport of the 
wave packet. Moreover, the direction of the transport is 
always perpendicular to and is thus controlled by the ex- 
ternal force. In view of the intensive investigations on the 
control of transport of ultracold atoms in optical lattices 
IL jl2i. il7i|. this fact may find use in future experiments. 
As for the experimental observation of the scenario con- 
sidered in this paper, we would say that the system can 
be readily realized with current experimental technology. 
The optical lattice can be constructed routinely 
the external force can be provided by the gravity 
tilting the two dimensional optical lattice from the hor- 
izontal plane, both the magnitude and the direction of 
the force can be adjusted at will. Furthermore, the great 
progress in in-situ observation of cold atoms in real space 
y\3\ guarantees that observation of the directed transport 
should also be achievable in the near future. Finally, we 
note that 1/Ji in Fig.[5l^a) corresponds to a time on the 
order of 0.6 ~ 0.7 ms for ^^Rb, while the life time of cold 
atoms in an optical lattice can be made well on the or- 
der of 1 s [ll|, [l2] ■ Thus this system is long-lived enough 
to allow for a significant drift which is in turn readily 
observable experimentally. 

An alternative approach for experimental realization 
is by the two dimensional photonic array. There, the 
propagation of light in the waveguides maps to the tight 
binding model perfectly, and the force can be effected by 
engineered refractive index variation Q. Actually, two 
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FIG. 1: (Color online) (a) Contour plot of the triangular optical lattice with the minima Vb = — 1.5i5r, with Er = /2m}? 
being the recoil energy. The basis vectors chosen are OA and OB, with \OA\ = \OB\ = 2A/3. The coordinates of sites A and 
B are (1, 0) and (0, 1), respectively, (b) Natural logarithm of the absolute values of the hoppings J(mi,m2) in units of Er- Note 
that the De symmetry of the lattice in (a) is respected, e.g., we have exactly J±(i,o) = ^±(o,i) = '^±(1,1) = ^^i- 
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FIG. 2: (Color online) (a) Time evolutions of the center-of-mass (mi) and (7712) of a wave packet in the optical lattice depicted 
in Fig.[T] The initial wave packet is a Gaussian tpm.{t = 0) = Aexp{ — {mf + m2)/a^ +ikQmi +ikQm2), with width a — 20 and 
wave vector fco = (fco, kg) — (0.05, 0.03). Evolution of the wave packet under three different forces are considered. They are (i) 
F/Ji = (0.5,-0.5), (ii) F/Ji = (0.7,-0.7), and (iii) F/Ji — (0.4,-0.8), respectively. In each case we have respectively, (i) 
(mi): blue 000 line, {m2): blue solid line; (ii) (mi): green + + + line, (7712): green dotted line; (iii) (mi): red * * * line, (m2): 
red dashed line, (b) Trajectories (mi) versus {m2) in the three cases. The solid, dotted, and dashed lines refer to the cases (i), 
(ii), and (iii), respectively. 



dimensional Bloch oscillation and Zener tunnelling have 
been observed in this system Therefore, it is likely 
that directed transport should also be observable in this 
system, which would show up as drift of the center of the 
pattern with the propagation length. 
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